Ordered Phase of the Dipolar Spin Ice under [110]-Magnetic Fields 
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We find that the true ground state of the dipolar spin ice system under [110]-magnetic fields is the 
"Q = X" structure, which is consistent with both experiments and Monte Carlo simulations. We 
then perform a Monte Carlo simulation to confirm that there exists a first order phase transition 
under the [110]-field. In particular this result indicates the existence of the first order phase transition 
to the "Q = X" phase in the field above 0.35 T for Dy2Ti207. We also show the magnetic field- 
temperature phase diagram to summarize the ordered states of this system. 
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Spin ice (SI) materials such as Dj^T^Oy and 
Ho2Ti207 have been attracting much interest as one of 
the ideal model materials of geometrically frustrated sys- 
tems 0, • In the SI materials, the magnetic moments 
occupy a three dimensional network of corner sharing 
tetrahedra with cubic symmetry called the pyrochlore 
lattice (FIG^I. Each magnetic moment is forced to point 
toward either of two tetrahedron centers by the strong 
crystal field along the local (111) axes 0,0. Then ferro- 
magnetic nearest neighbor interactions lead to a geomet- 
rical frustration among four spins in a tetrahedron and 
the local ground states in each tetrahedron are six-fold 
degenerated. These local ground states are determined 
by the rule that two spins point inward and two spins 
outward. This "2-in & 2-out" rule is called the ice rule 
from the analogy to the ice Ih [| . Under the ice rule, a 
number of low-lying states exist in the system and there 
remains a residual entropy almost same as the Pauling 
value HH0. 

It has been suggested that the large dipole interaction 
is responsible for the spin ice behavior in the SI materi- 
als Q ■ In the SI model with long range dipole interaction 
(the dipolar spin ice, DSI), it would be expected that the 
degeneracy of the ground state is removed by the long 
range nature of the dipole interaction and a long range 
order exists at low temperatures. Hertog et al. performed 
single spin flip Monte Carlo (MC) simulations of the DSI 
at zero field by using Ewald method and concluded that 
there is no long range order and there remains certain 
residual entropy close to the Pauling value Q . Melko et 
al, however, found the phase transition to the long range 
ordered q or d = (0, 0, lit /a) phase, where a is the size of 
conventional cubic unit cell, by MC simulations employ- 
ing the loop algorithm which improves the dynamics at 
low temperature. This q OT & order was identified as the 
true ordered ground state at zero field by the Fourier 
transformation of the interaction [Tol | . 

Recently, the SI materials under magnetic fields along 
the [110]-direction have been studied by means of sin- 
gle crystal experiments 0, Fennel et al. per- 
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FIG. 1: A schematic plot of the pyrochlore lattice and Q — X 
configuration. The magnetic ions are located at each vertex of 
the tetrahedra. The numbers in the figure indicate the Ising 
axis indices v — 0, 1, 2, 3. Under the [110]-fields, the system 
is divided into two sets of spin chains, a-chain (double line) 
and /3-chain (dotted line) . The [110]-field affects spins only 
on the a-chains. 



formed neutron diffraction experiments of Dy2Ti207 un- 
der the [110]-magnetic fields with fixed low tempera- 
ture ~ 70 mK, which shows magnetic Bragg peaks at 
the Q = points and diffuse scattering at the Q = X 
points From their scattering patterns, they found 

the coexistence of long-range ferromagnetic order on the 
field coupled spin chains (called a-chains by Hiroi et 
al. [Uj, FIGQJ and short-range antiferromagnetic order 
on the field independent spin chains (/3-chains). They 
suggested that the true ground state is the Q = X struc- 
ture as shown in FIG fTI [TL fTl| ■ although they could not 
find a clear evidence for such ordered state. A similar 
result was found in neutron diffraction experiments of 
Ho 2 Ti 2 7 by Harris et al Q. On the other hand, Hiroi et 
al. performed specific heat measurements of Dy2Ti207 in 
the [110]-magnetic fields. They found a relatively sharp 
peak at 1.1 K (below the broad peak associated with SI 
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freezing at low fields) in the field H > 0.4 T. They sug- 
gested that this peak comes from the freezing of /3-chains 
without long range order. They also argued the frustra- 
tion between /3-chains and threw doubt on the Q = X 
structure as the true ground state [l3( . In the theoretical 
point of view, Melko et al. performed MC simulations to 
observe a first order phase transition to the Q = X struc- 
ture ^3 • Thus the Q — X structure is still controversial. 

In the [110]-fields, there remains ground state degen- 
eracy on the /3-chains as far as only the nearest neighbor 
interaction is taken into account. We expect that this 
degeneracy is removed by the dipole interaction. In this 
paper, we give a proof that the Q — X structure is the 
true ground state configuration under the [110]-fields. We 
then perform single spin flip MC simulations to show that 
the relatively sharp specific heat peak originates from the 
first order phase transition to the Q = X phase and to 
present the temperature-field phase diagram. 

In the spin ice compound, each spin Si is forced to align 
along its easy axis rii due to the strong crystal field, so 
that it is expressed with an Ising variable Ui{— ±1) as 
Si = Uifii, where is one of the four distinct Ising axes 
n v \ v = 0, 1, 2, 3 corresponding to each site on the basic 
tetrahedron (FIG^I. Then the Hamiltonian of the Ising 
DSI system with N spins under an uniform external field 
H is given as 
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where the second and the last terms represent the near- 
est neighbor interactions (both exchange J nn and dipole 
Dnn) and the long range part of the dipole interactions, 
respectively. By using the size of unit cell a, the near- 
est neighbor distance r nn is written as r nn = a/>/8. 
In this paper we use material parameters of Dy2Ti207: 
J nn = —1.24 K, D nn = 2.35 K and p c sHB — 10.6/Xb = 
7.09 KT" 1 0. 

First, we consider the effect of magnetic field along the 
[110]-direction without the long range interaction term. 
In the [110]-fields, the pyrochlore lattice is divided into 
two sets of Ising spin chains named as a and /3-chains, 
as illustrated in FIGQ] [H E^. Spins on the v = 0, 3 
sites compose a-chains and ones on the v = 1, 2 sites 
/3-chains. The [110]-magnetic fields affect spins only on 
the a-chains and tend to align them ferromagnetically. 
Being free from the [110]-magnetic fields, spins on a (3- 
chain also align ferromagnetically due to the ice rule when 
the a-chains are ordered ferromagnetically. It should be 
noted that each ferromagnetic /3-chain has two possible 
states and the ground states are 2 Nf3 -fold degenerated, 
where Np ~ _/V 2 / 3 is the number of the /3-chains. At the 



ground state, the inter-chain configuration of a-chains is 
ferromagnetic (F) while the one of /3-chains is paramag- 
netic (P). Here we introduce abbreviated notations for 
inter-chain configurations of the a//3-chains for conve- 
nience. The inter-chain configuration of the a-chains is 
abbreviated as a-AF if they are ordered antiferromagnet- 
ically (AF). In this notation, the ground state configura- 
tion mentioned above is expressed as "a-F & /3-P" . Thus 
there is no phase transition to the long range ordered 
phase without the long range interaction. 

Now, we consider the effect of the long range interac- 
tion. At zero field, the long range interaction removes 
the macroscopic degeneracy of the ground states and 
the q or d = (0, 0, 2n/a) configurations become the true 
ground states [!J [Hj- I n terms of the inter-chain con- 
figuration of the a//3-chains, these ground states are ex- 
pressed as "a-AF & /3-AF". In order to show that the 
Q = X structure, which is expressed as "a-F & /3-AF" 
(see FIGQ, is the true ground state under [110]-magnetic 
fields, we only consider F and AF inter-chain alignments 
for the a//3-chains. By means of the Ewald summation 
techniques, we calculated the interaction energy of the 
four configurations for H = and obtained 
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Note that the two energy differences Eq~ f '''~ x> 



££- AF '' 3 ~ AF and E^~ F 



Eq~ F,13 ~ AF are both 



0.4122 K per spin. This result implies that there is 
no interaction between the a and /3-chains as far as F 
and AF inter-chain alignment are concerned. Indeed 
this statement is true since the internal fields acting on 
the /3-chains from the a-F/AF chains are canceled out 
each other due to the symmetry of the pyrochlore lattice 
structure. Since "a-AF & /3-AF" is the ground state 
at zero field, we can say that each set of ferromagnetic 
Ising spin chains tends to align anti-ferromagnetically. 
In the [110]-field, the energy of the "a-F & /3-AF" 
configuration comes down while the energy of the "a-AF 
& /3-AF" configuration does not change. So there is a 
ground state transition field H c , at which the following 
equation holds: 
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Using these values given in eq.J2J), we obtain the tran- 
sition point H c = 0.142 T. Thus we conclude that the 
"a-F & /3-AF" (Q = X) configuration is the true ground 
state in the field H > H c . 

Next we perform single spin flip MC simulations to 
show that there exists phase transition to the "a-F & 
/3-AF" phase. Our simulations were carried out on sys- 
tem size up to 5488 spins (cubic unit cell length L = 7) 
with a slow cooling process. The simulation lengths were 
~ 10 5 MC steps per spin at each temperature. In FIG. 
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Dy 2 Ti 2 7 ( L = 4, 1 024-spins ) 



Dy 2 Ti 2 7 ( L = 6, 3456-spins ) 
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FIG. 2: Temperature dependence of the specific heat for 
various magnetic fields along the [If 0]-direction. Each data 
set has been shifted for clarity. The dotted line indicates 
the zero level for each set. Four characteristic peaks at 
Tsi, Tap, T a , Tp are also indicated (see text). 



|2 we show our MC results for the temperature depen- 
dence of the specific heat for system size L — 4 under 
various magnetic fields H . At zero field a broad peak 
corresponding to the freezing into the SI state is found 
around Tgi ~ 1.2 K, in good agreement with the previous 
results 8j. As the magnetic field increases, the peak at 
T S i grows up to H = 0.30 T, and then at H = 0.35 T a 
small anomalous peak appears at T a p on the low temper- 
ature side of Tgi, and this anomalous peak moves to the 
high temperature side. At H — 0.50T, the peak at T a p 
merges with the broad peak at T$i and there appears a 
very sharp peak, and then this peak is broken into two 
peaks at T a and Tp. The peak at T a is broad and shifts to 
the high temperature side with increasing field. On the 
other hand, the peak at Tp is sharp and almost indepen- 
dent of the field. From their field dependences, the peaks 
at T a and at Tp must come from the a and /3-chains, re- 
spectively. To extract the contribution of the /3-chains, 
we performed a MC simulation of the system only with 
spins on the /3-chains. This system equals to the DSI at 
the high-field limit since the spins on the /3-chains feel no 
internal field from the a-F chains. The temperature de- 
pendence is also shown in the FIG.|21and we see a sharp 
peak near Tp ~ 1.3 K. This sharp peak corresponds to 




FIG. 3: Temperature dependence of the /3-AF order param- 
eter for [110]-magnetic fields H = 0.35 T, 0.50 T and 1.00 T . 
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FIG. 4: Finite size scaling of the peak height of the specific 
heat as a function of system size N = 16L 3 for L=3,4,5,6,7. 



the peak at Tp in large but finite magnetic fields as ex- 
pected. These temperature and field dependences of the 
specific heat are also seen in the experiments by Hiroi et 
al\jj^. 

In order to consider what happens at each peak more 
precisely, we observe order parameters corresponding to 
the inter-chain orders. For example, the order parameter 
for the /3-AF order mp-AF is the q rd = (0, 0, 2-k/o) stag- 
gered magnetization of the /3-chains. In FIG|2lwe show 
the temperature dependence of mp.^F for system size 
L = 6. We can see the /3-AF order below T a p and Tp, 
and discontinuity of the order parameter mp.^F at these 
temperatures. We also find from the observation of m a .F 
that the field sensitive broad peak at T a comes from the 
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FIG. 5: Magnetic field-temperature phase diagram for H || 
[110]. The open squares, closed squares, open circles and 
closed circles indicate the characteristic peaks Tsi, T a p, T a 
and T/3, respectively. The solid line shows the first order phase 
transition and the dotted the crossover. 



a-F freezing. Therefore the low temperature anomalous 
peak at T a p and the field independent peak at Tp are due 
to the first order phase transition to the "a-F & /3-AF" 
ordered phase. We also carried out a finite size scaling of 
the peak height C pea k for the magnetic field H = l.OOT 
and for the /3-chain system (FIG0J. The scaling data 
shows that C pea k is proportional to the system size I? 
for large systems (L > 5). These results prove that the 
phase transition is indeed of the first order. 

Finally we summarize our results from the ground state 
analysis and the MC simulation as a phase diagram in 
FIGEJ The data points in the diagram indicate the spe- 



cific heat peak temperatures Tgi, T a p, T a , Tp in the MC 
results with single-spin-flip slow cooling process (FIGJ2J- 
The solid line in the diagram corresponds to the first or- 
der phase transition. Note that the straight lines shown 
in the low temperature and in the weak field region are 
expectation ones from the ground state transition field 
H c = 0.142 T and the zero field MC results @. In the 
cooling process at high field region, spins on the a-chains 
are almost fixed by the [110]-magnetic field below T a , 
while the spins on the /3-chains still take random con- 
figurations ("a-F" phase). Then at Tp the spins on the 
/3-chains form the /3-AF order and the phase transition 
to the "a-F & /3-AF" configuration takes place due to 
the ice rule and the long range dipole interaction. In the 
weak field region H < 0.30 T, the phase transition has 
not been observed in our single spin flip MC simulation. 

In conclusion, we have studied the effect of the dipole 
interaction on the DSI system under the [110]-magnetic 
fields through the ground state analysis and the single 
spin flip MC study. We have found that the "a-F & 
/3-AF" (Q — X) configuration is the true ground state in 
the strong field region H > H c = 0.142 T and the first 
order phase transition from the "a-F" phase to the "a-F 
& /3-AF" phase occurs. We have presented the [110]-field- 
temperature phase diagram. Melko et al. applied the loop 
algorithm and the multicanonical MC method in order to 
avoid the freezing at low temperature in the first order 
phase transition [3, Il2| . Such techniques will clarify the 
detail of phase boundary in the low temperature and the 
weak field region, and are now under progress. 

We are grateful to Z. Hiroi, K. Matsuhira, M. Ogata, 
and T. Sakakibara for useful comments and discussions. 
This research was partially supported by the Ministry 
of Education, Science, Sports and Culture, Grant-in- Aid 
for JSPS Fellows. 



[1] M. J. Harris, S. T. Bramwell, D. F. McMorrow, T. Zeiske 
and K. W. Godfrey: Phys. Rev. Lett. 79 (1997) 2554. 

[2] S. T. Bramwell and M. J. Gingras: Science 294 (2001) 
1495. 

[3] R. Siddharthan, B. S. Shastry, A. P. Ramirez, A. 

Hayashi, R. J. Cava and S. Rosenkranz, Phys. Rev. Lett. 

83 (1999) 1854. 
[4] S. T. Bramwell, M. N. Field, M. J. Harris and I. P. 

Parkin, J. Phys. Condens. Matter 12 (2000) 483. 
[5] A. P. Ramirez, A. Hayashi, R. J. Cava, R. Siddharthan 

and B. S. Shastry: Nature 399 (1999) 333. 
[6] L. Pauling : J. Am. Chem. Soc. 57 (1935) 2680. 
[7] S. Yoshida, K. Nemoto and K. Wada: J. Phys. Soc. Jpn 

71 (2002) 948. 



[8] B. C. den Hertog and M. J. P. Gingras: Phys. Rev. Lett 

84 (2000) 3430. 
[9] R. G. Melko, B. C. den Hertog and M. J. P. Gingras: 

Phys. Rev. Lett. 87 (2001) 067203. 
[10] M. J. P. Gingras and B. C. den Hertog: Can. J. Phys. 79 

(2001) 1339. 

[11] T. Fennell, O. A. Petrenko, G. Balakrishnan, S. T. 

Bramwell, J. D. M. Champion, B. Fak, M. J. Harris and 

D. M. Paul: Appl. Phys. A 74 (2002) S889. 
[12] R. G. Melko, M. Enjalran, B. C. den Hertog and M. J. 

P. Gingras: cond-mat / 0308282. 
[13] Z. Hiroi, K. Matsuhira and M. Ogata: J. Phys. Soc. Jpn 

72 (2003) 3045. 



